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We consider the stochastic background of gravitational waves produced during the radiation- 
dominated hot big bang as a constraint on the primordial density perturbation on comoving length 
scales much smaller than those directly probed by the cosmic microwave background or large-scale 
structure. We place weak upper bounds on the primordial density perturbation from current data. 
Future detectors such as BBO and DECIGO will place much stronger constraints on the primordial 
density perturbation on small scales. 



I. INTRODUCTION 



o 

Recent cosmic microwave background (CMB) experiments [1, 2] and high-redshift galaxy surveys [3] are able to 
probe the primordial density perturbation on cosmological scales, ~ 10 — 1000 Mpc. This yields a precise measure 
of the primordial density contrast, A^(fcc , A/s) = 2.5 x 10~ 9 for kcMB — 0.002 Mpc -1 [1]. For much smaller scales, 
7-H ' free-streaming of relativistic particles in the very early universe and Silk damping due to photon diffusion erases the 
primordial density perturbation on comoving scales much less than 10 Mpc [4]. The present distribution of matter on 

O smaller scales is the result of the subsequent non-linear evolution of the matter density and we are unable to directly 
relate observational data to the primordial distribution of matter. Therefore for modes which are much smaller than 
u: 10 Mpc the primordial density perturbation remains undetermined. 

One might assume that the primordial density perturbation on small scales has the same value as on larger scales. 
O ,, According to WMAP 5-year data [1, 2] the spectral index on CMB scales is smaller than one, n s = 0. 96^ooi3) 5 
implying that the density contrast is slightly smaller on smaller scales, but it is a huge extrapolation to assume that 
this is the actual scale-dependence all the way from CMB scales down to comoving scales such as < 10 pc, for 
instance, which are smaller than the comoving Hubble scale at the epoch of big bang nucleosynthesis, T ~ MeV. The 
| only probes of such small scales at early times are expected to be gravitational relics. 

One process which allows us to put an upper bound on density perturbations is the formation of primordial black 
| holes (PBHs) [5]. PBHs are produced when density fluctuations with a large amplitude (Sp/p ~ 0.01 — 0.1) enter the 
, horizon. Their typical mass is given by the horizon mass when they were produced [6-9]. 

-t ; . . . 

^ | , which can be related to the temperature if they formed in the radiation dominated era 

SI M PBHs *10™^f)\ (2) 

• The fact that PBHs have not been observed to date limits the initial mass fraction going into black holes, which can 

be used to constrain the primordial density perturbation on the corresponding scales [6]. For instance PBHs radiate 
Hawking radiation and thus for Mpbhs < 10 15 g evaporations are limited because of constraints from the standard 
big bang nucleosynthesis (BBN), which is sensitive to the baryon-photon ratio at the time, and also because of bounds 
on the expected gamma-ray background. For heavier PBHs (Mpbhs > 10 15 g), similar constraints come with the fact 
that the present-day density of PBHs cannot exceed the upper limit on the cold dark matter density [10]. 

Gravitational waves also give us a window onto density perturbations in the very early universe. They can be 
generated by violent events in the early universe such as bubble collisions [11], cosmic strings [12], preheating after 
inflation [13-16], or parametric decay of supersymmetric condensates [17, 18]. All of these mechanisms are highly 
model-dependent. Even gravitational waves from produced from vacuum fluctuations of the metric during inflation 
[19-21], though a generic prediction of inflation, have an amplitude which is dependent upon the energy scale of 
inflation. 

In this paper we consider the bounds placed on the the primordial density perturbation from the generation of a 
stochastic background of gravitational waves in the very early universe. In the standard radiation-dominated hot big 
bang, first-order density perturbations inevitably generate gravitational waves at second (and higher) order [22-32] . 
Like the initial mass fraction of PBHs, the power spectrum of induced gravitational waves is determined by the 
primordial density perturbation and thus can place bounds on the amplitude of scalar perturbations. Indeed Saito 
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and Yokoyama [34] have recently pointed out that current gravitational wave bounds are sufficient to rule out PBHs 
as possible candidates for intermediate mass black holes. 

If we could determine the amplitude and frequency of these induced gravitational waves then we would be able 
to determine the primordial density perturbation when the corresponding scales crossed the Hubble scale during the 
early hot big bang. Throughout this paper we will use the following formula to relate the frequency of gravitational 
waves at the present time to the temperature at Hubble-crossing in the early radiation-dominated era [33] : 

.-i-UxIO-V (3) 

where g* is the effective number of degrees of freedom. We are able to place upper bounds on the primordial density 
perturbation from BBN and cosmic microwave background (CMB) constraints and as well as current LIGO, VIRGO 
and pulsar timing data. If we do not detect gravitational waves with future pulsar timing [35, 36] and future detectors 
like Advanced LIGO [37], Advanced VIRGO [38], LISA [40], BBO [42] and DECIGO [43], the upper bounds on the 
primordial density perturbation will become significantly tighter. 

We emphasize that our bounds come from adopting the standard, minimal cosmological model of adiabatic density 
perturbations, in their growing mode, in a radiation-dominated early universe from ultra-high energies (~ 10 16 GeV) 
until matter-domination when T < eV. The quantitative constraints will be altered if one adopts non-standard 
cosmological evolution [44, 45] such as an early matter-era (e.g., temporary domination of the energy density by 
massive, non-relativistic particles) or a stiff- fluid-dominated era (e.g., domination by the kinetic energy of a coherent, 
fast-rolling scalar field). On general grounds one expects an early era dominated by fluid "softer" than radiation, 
P/p < 1/3, to dilute the fractional density of gravitational waves whose wavelength is smaller than the comoving 
Hubble scale (and thus behave like radiation) , while the fractional density of sub-horizon gravitational waves grows 
relative to matter which is stiffer than radiation [47]. On the other hand the evolution of density perturbations 
which give rise to gravitational waves is also altered, see for example Ref. [33]. Non-adiabatic modes in a multi- 
component system can lead to large-scale adiabatic density perturbations by the time of last-scattering but arise 
from initial isocurvature perturbations, so are not necessarily constrained by our analysis. One such example is the 
curvaton scenario, where the gravitational waves may be enhanced with respect to the adiabatic case if the curvaton 
is subdominant when it decays [48]. Such models need to be considered on a case by case basis. 

This paper is organized as follows: in section 2 we introduce the basic equations used to determine the induced 
gravitational wave background and define the effective energy density of second-order gravitational waves. In section 3 
we quantify the constraints placed on the primordial density perturbation by a variety of experiments. We present 
our conclusions in section 4. 



II. SECOND-ORDER GRAVITATIONAL WAVES 



In this section we will briefly review the generation of induced gravitational waves. Details of the calculations have 
been described previously [30, 31]. 

The perturbed metric in the longitudinal gauge is 

ds 2 = a 2 (7?)[-(l + 2<S>)d V 2 + [(1 - 2*)% + 2F {iii) + h iJ ]dx i dx i ] (4) 

where <j> and <]/ are scalar metric perturbation, Fi is a transverse vector perturbation and hij is a transverse and 
trace-free tensor perturbation. The scalar metric perturbations, $ and are supported by density perturbations, 
and in the absence of anisotropic stress we require $ = * [49] . We will find it convenient to use the Fourier transform 

$ ( X ) = 7-TT /d 3 k<iW kx , (5) 



(27T) 

where, for an isotropic distribution, the power spectrum is given by 

2tt 2 

($k$k') = -p-S 3 (k + kOP(/c) (6) 

On large scales (much larger than the Hubble scale) the power spectrum of the primordial scalar perturbation is 
commonly approximated by a power law 



V{k) = -A£ ( - ) (7) 

. ft* 
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where the numerical factor 4/9 comes from the relation between scalar metric perturbation in the longitudinal and 
comoving gauges on large scales in a radiation-dominated era [49]. 

In the second-order perturbed Einstein field equations we see the effect of first-order perturbations as a source term 
(Sij) for second-order tensor perturbations. After putting all first-order perturbation terms to the right-hand side of 
the Einstein field equation, we have 

+ ZHh'y + k 2 h tl = Sl T , (8) 

where S^ T indicates the transverse-tracefree part of the source term. If we neglect first order tensor and vector 
perturbations in comparison with first order density perturbations, the right hand side of this equation is the trans- 
verse and trace-free part quadratic in first-order scalar perturbations. This behaves like a source term for induced 
gravitational waves [30, 31] 

S,^ = 2$didj$ - 2*<9 i <9 7 $ + 4*<9 i d J * + d^dj® - d^djV - d'Vdj® + 38^8^ 
2c 2 



3wU 



[3H{H^-^') + V 2 ^}d i 8 J (^-^) (9) 



where w = P/p is the equation of state and c 2 = P'/p' is the adiabatic sound speed. 

These equations are written in the real space but in order to derive the power spectrum of gravitational waves we 
need to transform to Fourier space [30] 

/d 3 k 
— T e ,k - x [/i k (7 7 )e ij -(k) + fckey(k)] , (10) 

where e y (k) and ejj(k) are the polarization tensors. The two polarization tensors etj(k) and eij(k) can be given in 
terms of the orthonormal basis 

e»j(k) = ^[e»(k)ej(k) -ei(k)ej(k)] 

ey(k) = -J=Mk)e J (k) + e,(k) ej (k)], (11) 

where e and e are unit vectors orthogonal to one another and k: 

ei k l = e.k' 1 = e.e 1 = (12) 
The gravitational waves have a power spectrum in Fourier space 

1 2-7T 2 

(h^v)hv(ri)) - --^(k + kOTMM), (13) 

The effective density of a stochastic background of gravitational waves, on scales much smaller than the Hubble scale, 
is given by [21] 

p ™ - aib^ - 32^ S d ^ h) Vh ^) • (14) 

The fraction of the critical energy density in gravitational waves per logarithmic range of wavenumber k in the 
radiation era is thus 

^W(M) = ^(£) Ph(k,v). (15) 

After the radiation-dominated era, the density of gravitational waves on sub-Hubble scales then redshifts exactly as 
any non-interacting relativistic particle species and in the present day we have 

n c^.oW = % 2 (^) n(k, V ), (16) 

where we neglect additional numerical factors due to the detailed thermal history, such as the heating of photons by 
the annihilation of other relativistic particle species [45, 50]. The present density of photons is J7 7i o — 4.8 x 10 -5 
where here, and throughout this paper, we take H ~ 72 km s _1 Mpc _1 for the present value of the Hubble rate. 
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III. CONSTRAINTS ON PRIMORDIAL DENSITY PERTURBATIONS 



In the standard cosmological scenario, second-order gravitational waves are generated during the radiation- 
dominated era after inflation. Non-linear interactions can in principle lead to density perturbations integrated over 
a range of scales contributing to the gravitational wave amplitude on a given wavenumber, k, but in practice the 
second-order gravitational waves are primarily produced when first-order density perturbations on the similar on 
same scale, <~ k, come inside the Hubble scale during the radiation era [30]. 

Assuming a power-law spectrum for the primordial density perturbation, Eq. (7), the energy density of second-order 
gravitational waves, relative to the critical density at the present time, which were produced during the radiation- 
dominated era [y > 10~ 15 Hz), can be written as 

n gw . (k) = F rad n 1 . A i n (k). (17) 

where, for modes which are well inside the horizon at the end of the radiation-dominated era ikr\ eq » 1), we have [30] 

8 /216 2 \ 

F rad = 3 ( ) 8 ' 3 X 10 " 3 /ns ( 18 ) 

and f ns is weakly-dependent on the spectral tilt. f ns « 1 if the density perturbation is scale-invariant [30], but 
becomes slightly smaller than one for a red spectrum (e.g, f ns w 0.97 for n s = 0.9) and bigger than one for a blue 
spectrum (e.g, f ns w 1.05 for n s = 1.1). 

It is also possible to consider the spectrum of gravitational waves generated by density perturbations with a sharply 
peaked power spectrum [30, 34]. Considering a delta-function power spectrum, P(k) — (4/9)A^(fc p )<5(ln(fc/fc p )), the 
resulting gravitational wave spectrum is described by a sharply rising spectrum for k < k p [34] 

n flW ,o(fe) = 29fi 7i oA^(fc p )^ , (19) 

with an abrupt cut-off for k > 2k p . 

In the following numerical estimates we take F rad w 30 in Eq. (17) corresponding to an approximately scale- 
invariant spectrum of scalar perturbations, n s w 1. This is expected to be a conservative lower bound on F rac i for 
the blue spectra with n s > 1 required to produce a detectable background of gravitational waves on scales much 
smaller than the CMB scale. In the rest of this section we show, how Eq. (17) enables us to use constraints on the 
stochastic background of gravitational waves generated during the radiation era, Slg W ,o(k), to place upper bounds 
on the primordial density perturbation on the corresponding scales, A^(fc). In addition Eq. (19) indicates how 
observational constraints on {l gw fi{k) at a given wavenumber k also places a weaker bound on the primordial density 
perturbation, A^(fc p ) oc (k p /k), at higher wavenumbers, k p > k. Our results are presented graphically in Figure . 



A. Cosmological density constraints 

1. Constraint from BBN 

If the energy density carried by gravitational waves at the time of primordial big bang nucleosynthesis (BBN) were 
large, the abundances of the light nuclei produced would be altered with respect to the predictions of standard BBN. 
Hence, BBN can be used to constrain the total energy carried by gravitational waves at the time of nucleosynthesis 
(T~ 1 MeV) [51]. 

Primordial abundances of the light elements, usually quoted as a bound on the effective number of relativistic 
species at the time of BBN, gives the 95% c.l. upper bound on a primordial gravitational wave background [37] 

JW < 1.5 x 10~ 5 (20) 

Substituting this bound into Eq. (17) gives 

A 2 n < 0.1 

This denotes the upper bound on the primordial density perturbation on the Hubble scale at the time when the 
gravitational waves are generated. 



Fraa 
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FIG. 1: Constraints on the primordial density perturbation, A^, obtained from gravitational waves produced during the 
radiation era, using Eq. (17). Black lines denote current constraints from gravitational waves detectors and BBN. Green lines 
denote constraints expected from future gravitational waves detectors. Diagonal lines show the bounds on A^j obtained using 
Eq. (19) for wavenumbers k p > k. WMAP gives a direct measurement (shown in red) of the primordial density perturbation 
on very low frequencies. 



Although Eq. (21) is only a weak limit on the primordial density perturbation, it is applies across a wide range of 
length scales. It applies on scales which are smaller than the Hubble scale at the epoch of BBN, which from Eq. (3) 
corresponds to frequencies v > 10~ 10 Hz today, and scales which are larger than the Hubble scale at the start of the 
radiation-dominated era. This is mo del- dependent, but in an inflationary cosmology this would be the Hubble scale 
at the end of reheating after inflation, which could be as large as v ~ 10 8 Hz for T ~ 10 16 GeV. If inflation occurs at 
lower energy scales the reheating temperature could be much lower. 

The relationship between the primordial density perturbation on the CMB scale, /\\{hcMB)i and on an arbitrary 
scale, A^.(fc), can be written as 

,2 



A^(fc) = A^(k CMB ) x (22) 

where n s describes the effective tilt between the scale k and CMB scales, where is directly observed for 
A^(A;cms) = 2.5 x 10~ 9 for kcMB = 0.002 Mpc -1 [1]. In terms of frequency we have 



log 



' ■ L 



log ( —*—) 



(23) 



where vcmb ~ 10~ 18 Hz. Therefore (21) can be interpreted as a constraint on the value of the effective spectral 
index, n a , across a wide range of scales. 

According to WMAP5 data [1], the spectrum of primordial density perturbations is red (n s less than unity) on 
CMB scales but, except for the constraints from primordial black holes (PBHs), we have no restrictions on the value 
of n s on scales much smaller than 1 Mpc. We obtain the tightest constraint on n s from a bound such as Eq. (21) 
applied to the smallest possible wavelength. For v ~ 10 8 Hz from (21) and (23) we find 



^<1.29-^log 10 (^). (24) 

Equation (24) corresponds to the maximum possible reheat temperature after inflation. However if we consider a 
smaller reheat temperature, for instance T ~ 10 s GeV corresponding to v ~ 1 Hz, we get n s < 1.44— 4^ log 10 
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2. Constraint from CMB 

A very similar bound on the effective energy density of primordial gravitational waves can be obtained around 
the time of last scattering of the cosmic microwave background. Again a limit on the number of massless neutrino 
species [52] can be translated into a bound on the gravitational wave background [53]. Unlike the BBN constraint, 
however, the CMB constraint depends upon the nature of inhomogeneous perturbations about the average density. 
For a gravitational wave background produced from a Gaussian random field of primordial density perturbations on 
small scales, we expect the effective energy density on long wavelengths (on scales of order 100 Mpc) to be independent 
of the density perturbations on this scale. Thus long wavelength perturbations of the gravitational wave background 
are non-adiabatic and Smith et al [53] give a 95% c.l. bound 

n gwfi < 1.3 x 10- 5 , (25) 

for a "homogeneous" gravitational wave background. This is marginally stronger than the BBN constraint (20). It 
gives effectively the same constraint on the primordial density perturbation (21), and the effective spectral index 
(24), but extends to longer wavelengths <~ 10~ 15 Hz, corresponding to scales inside the Hubble scale at the time of 
last-scattering. 

Future data from CMB experiments such as Planck and the proposed CMBPol mission are expected to improve 
the CMB bound. For Planck the expected bound corresponds to Q gw fl < 2.7 x 10~ 6 [53] which would bound 

_ i_ 

A^<0.04(^)"\ (26) 

B. Constraints from ground-based detectors 

1. Current LIGO/VIRGO 

We can obtain a tighter constraint on the primordial density perturbation on scales probed by direct detectors, 
such as the Laser Interferometer Gravitational Wave Observatory (LIGO) [37] and gravitational wave detector at the 
European gravitational observatory (VIRGO) [38]. LIGO's maximum sensitivity is around a frequency, v = 100 Hz. 
The latest results from the LIGO S5 science run give a bound on the energy density of gravitational waves on this 
scale [39] 

fi ff «,,o < 6.9 x 10~ 6 . (27) 
Hence from Eq. (17) the constraint on the density perturbation on the LIGO/VIRGO scale is 

A^<0.07(^)"\ (28) 

This is a slightly tighter bound than that currently obtained from BBN and the CMB, Eq. (21), however unlike the 
BBN and CMB bound it only applies to LIGO/VIRGO scales. The corresponding constraint on n s on this scale 
comes from Eq. (23) 

Ti a <1.37-llog 10 (^) . (29) 
2. Advanced LIGO/VIRGO 

Advanced LIGO/VIRGO will give us an improved constraint on a stochastic background of gravitational waves 
on the same scales [38, 54]. The smallest density of gravitational waves which could be detected by Advanced 
LIGO/VIRGO is 10 3 times smaller than current LIGO/VIRGO bounds. Considering the smallest detectable energy 
density f2 ffu > : o < 10~ 9 in Eq. (17) returns 

_ i 

A^<8xHT 4 (^) 2 (30) 
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Equation (23) gives the expected constraint on n s from Advanced LIGO/VIRGO (taking v = 100 Hz) 

*<i-n-Sj>*.(^)- (3D 

C. Constraints from LISA 

The Laser Interferometer Space Antenna (LISA) is the first gravitational wave detector planned in space and is 
the most sensitive detector currently planned. Assuming LISA's instrumental noise is well-behaved [55], it could 
detect a stochastic background of gravitational waves at a level Q gw fi ~ 10~ n at frequencies vlisa ~ 10~ 3 Hz 
[40, 41]. However the sensitivity of LISA leads to many potential overlapping sources and hence the problem of source 
confusion. In particular the astrophysical background from unresolved extra-galactic white-dwarf binaries is expected 
to limit LISA's ability to distinguish a primordial gravitational wave background to [55, 56] 

tt gw ,o < 1(T 10 . (32) 
The corresponding upper bound on the primordial density perturbation on LISA scales comes from (17): 

A^<3xl0- 4 ^) 2 (33) 
The constraint on n s on LISA scales comes from Eq. (23) 

W.<1.34-llo gM (%i). ,34) 

This is a slightly weaker bound on the effective spectral index compared with Advanced LIGO, as LISA is sensitive 
on length scales much larger than LIGO scales. 

D. Constraints from BBO/DECIGO 

The Big Bang Observer (BBO) [42] and the DECi-hertz Interferometer Gravitational wave Observatory (DECIGO) 
[43] are ambitious proposals for future space-based observatories to detect cosmological gravitational waves. They 
should be able to detect a stochastic background of gravitational waves down to an effective energy density fl gw .o ~ 
10~ 16 at vb bo ~ 1 Hz. This waveband is chosen to avoid the confusion noise due to white dwarf binary mergers 
which cuts off above 0.2 Hz. The designs of BBO and DECIGO are based on the requirement to identify and remove 
the remaining foregrounds from neutron star and black hole binaries [57]. 

If induced gravitational waves during the radiation era are not detected with BBO/DECIGO, then we will be able 
to place a tight constraint on the primordial density perturbation and hence n s on this scale (1 Hz). From (17) and 
(23) we obtain 

< 3 x 10~ 7 

and hence 

E. Constraints from pulsar timings 

Analysis of pulse data from pulsars shows that they are very stable clocks. Measurement of timing residuals, 
which is the difference between the observed time of arrival and predicted time of arrival, can in principle be used 
to directly detect gravitational waves passing between the pulsar and the observer [35, 36]. The data from current 
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observations of an array of pulsars places an upper bound on the stochastic background of gravitational waves, with 
periods comparable to the total observation time span. This is typically 1-10 years, and corresponds to 10~ 8 -10~ 9 
Hz. For v = 10~ 8 Hz, the constraint on the present density of gravitational waves is [36] 

O ffm ,o<4xl0- 8 . (37) 

Substituting (37) in (17), gives us the current constraints on the primordial density perturbation 

_ i 

A^<5xl0- 3 (^f\ (38) 

The constraint on n s comes from (17) 

To 10 *" {¥)■ (39) 

Saito and Yokoyama [34] have recently used similar constraints, on the induced gravitational wave background 
from pulsar timing arrays, to rule out the large amplitude of primordial density perturbations required to produce 
any significant number of primordial black holes in the intermediate mass range, 4 x 10 2 M Q < Mpbh < 5 x 10 3 M q , 
corresponding to 8 x 10 35 g < Mpbh < 10 37 g, which from Eq. (2) would have formed at temperatures T w 3—10 MeV. 

Future pulsar timing will give a better constraint. If gravitational waves are not detected, the upper limit, based 
on timing 20 pulsars over 5 years, would be Q gw fl < 10~ 10 [36]. From (17), the future constraint on the primordial 
density perturbation in five years time for v = 10~ 8 Hz would be 

/\\ < 3 x 10~ 4 

The future constraint on n s would then be 

".<1.50-|;log, (^). (41) 
IV. CONCLUSION 

Despite remarkable recent progress in astronomical observations mapping density perturbations on large scales (10- 
1000 Mpc) in our Universe, we know little about the primordial distribution of matter on much smaller scales. This 
is due to Silk damping and free-streaming of relativistic particles in the early universe, and subsequent non-linear 
evolution of matter perturbations at much later times. The only constraints on these scales come from gravitational 
relics of the very early universe. Previous work has focussed on the possible formation of primordial black holes from 
large over-densities. 

In this paper we have shown how limits on a stochastic background of gravitational waves can be used to place 
limits on density perturbations in the early radiation-dominated era of the standard hot big bang cosmology. 

BBN and CMB limits on a primordial gravitational wave background places only a weak constraint on the amplitude 
of primordial density perturbations, < 0.1, but this applies across a wide range of frequencies, from 10~ 15 Hz 
to frequencies as high as 10 8 Hz, depending on the maximum temperature at the start of the radiation-dominated 
era. By contrast, gravitational wave detectors such as LIGO and VIRGO place slightly tighter bounds, currently 
< 0.07, but only over a narrower range determined by the frequency response of the detector. 

Future gravitational wave experiments offer the prospect of much tighter bounds on, or a detection of, a stochastic 
gravitational wave background and hence the primordial density perturbation on small scales. A space-based experi- 
ment such as LISA could detect gravitational waves produced by density perturbations A^ ~ 10~ 4 , and future data 
from pulsar timing arrays could have similar sensitivity. The most ambitious current proposed gravitational wave 
observatories including BBO and DECIGO offer the prospect of detecting gravitational waves as small as A|j ~ 10~ 7 . 

If gravitational wave background is not detected by these experiments it would imply that the primordial power 
spectrum remains close to scale-invariant, or decreases in power on small scales, n s < 1.29, which provides a valuable 
new constraint on models for the origin of structure. Nonetheless it remains a challenge to design an experiment that 
could detect gravitational waves produced by primordial density perturbations of the same power, A^ ~ 10~ 9 , as 
seen on the largest scales in the universe today. 



/ Frad\ 
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We have assumed the simplest expansion history of the universe, being radiation dominated from very early times, 
corresponding to temperatures as high as 10 16 GeV. If the early universe has a more complicated history, the con- 
straints may be altered. A period of inflation is expected to dilute pre-existing gravitational waves on sub-Hubblc 
scales, without generating a significant additional background [58]. On the other hand although an early matter- 
dominated era before BBN, such as the reheating or preheating after inflation, would also dilute any gravitational 
waves that had already been generated, it could itself produce significant tensor metric perturbations on scales that 
re-enter the Hubble scale during an early matter-dominated era [33]. We leave a more detailed investigation of 
constraints on the primordial density perturbation in more general cosmological scenarios to future work. 
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